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Doelstelling 

Het  doel  van  dit  project  was  het  vergroten  van  de  inzetbaarheid  van  SAR 
(Synthetische  Apertuur  Radar),  en  PHARUS  (vliegtuig  SAR)  in  het  bijzonder,  ten 
behoeve  van  de  detectie  en  herkenning  van  niet  bewegende  militaire  doelen.  Het 
accent  lag  hierbij  op  detectie-  en  herkenningsalgoritmen  die  gebruik  maken  van 
polarimetrische  SAR  data. 

Omschrijving  van  de  werkzaamheden 

In  de  eerste  fase  van  dit  project  is  een  cursus  gevolgd  op  het  gebied  van  detectie  en 
herkenning  en  is  een  grote  literatuurstudie  uitgevoerd  naar  detectie-  en 
herkenningsalgoritmen  voor  polarimetrische  SAR  data.  De  meest  belovende 
algoritmen  zijn  vervolgens  geimplementeerd  en  getest  op  PHARUS  beelden  van 
Swynnerton  (UK),  opgenomen  tijdens  de  NAVO  Panel  3  RSG-20  MIMEX 
proefnemingen  in  oktober  1996.  De  kwaliteit  van  deze  beelden  is  grondig 
onderzocht  ten  behoeve  van  de  evaluatie  van  de  resultaten. 

Resultaten 

SAR  data  bevat  meer  mis  dan  optische  of  infrarood  beelden.  Om  het  aantal  false- 
alarms  bij  detectie  en  herkenning  te  beperken,  wordt  deze  mis  vaak  uitgefilterd. 
Diverse  polarimetrische  filter  algoritmen  zijn  onderzocht.  Het  beste  algoritme 
bleek  het  Polarimetrisch  Whitening  Filter  (PWF).  Het  detecteren  van  militaire 
voertuigen  met  PHARUS  bleek  vervolgens  zeer  succesvol.  Van  de  onderzochte 
detectoren  bleek  de  Order-Statistics  Constant  False-Alarm  Rate  (OS-CFAR) 
detector  het  meest  robuust.  Het  herkennen  van  militaire  voertuigen  op  basis  van  de 
vorm  en  de  polarimetrische  eigenschappen  bleek  niet  mogelijk.  De  belangrijkste 
oorzaak  hiervan  was  de  te  lage  resolutie  van  PHARUS  die  hier  3.5  m  x  5.0  m 
bedroeg.  In  het  vervolg  op  dit  onderzoek  [de  Vries  1998]  is  gekeken  naar  PWF 
gefilterde  data  met  een  hogere  resolutie  dan  PHARUS  op  dit  moment.  Het  strekt 
tot  de  aanbevelingen  om  ook  te  kijken  naar  de  polarimetrische  eigenschappen  van 
hoge  resolutie  data.  Detectie  en  herkenning  kan  het  beste  worden  uitgevoerd  in 
combinatie  met  andere  sensoren  als  thermisch-infrarood  (TIR),  laser-radar  en 
multispectrale  scanners.  Het  voordeel  van  SAR  en  laser-radar  is  dat  zij  ook  in  staat 
zijn  hoogte  of  diepte  weer  te  geven. 
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1. 


Inleiding 


Voor  het  opsporen  van  doelen  kunnen  verschillende  sensoren  worden  gebruikt. 
SAR  (Synthetische  Apertuur  Radar)  is  aan  te  bevelen  omdat  deze  ongeacht  het 
weer  en  ongeacht  of  het  dag  of  nacht  is,  luchtopnamen  kan  maken.  TNO-FEL  heeft 
een  eigen  SAR  systeem  PHARUS  (Phased  Array  Universal  SAR)  [Hoogeboom 
1994]  ontwikkeld,  in  samenwerking  met  NLR  en  TU  Delft.  Het  bijzondere  aan  dit 
systeem  is  dat  het  volledig  polarimetrisch  is  en  daardoor  meer  informatie  verza- 
melt  over  de  aard  van  de  belichte  objecten.  Daamaast  is  PHARUS  uitgerust  met 
een  phased  array  antenne,  zodat  het  mogelijk  is  de  richting  van  de  antenne  bundel 
te  varieren  tijdens  de  vlucht,  wat  de  inzetbaarheid  vergroot. 

Dit  rapport  is  een  verslag  van  onderzoek  dat  is  uitgevoerd  naar  methoden  om 
doelen  op  te  sporen  met  SAR,  en  in  het  bijzonder  met  PHARUS.  Een  doel  kan 
hierbij  ieder  persoon,  object  of  verschijnsel  zijn  waarin  men  is  gei'nteresseerd.  Bij 
het  opsporen  van  doelen  kunnen  de  volgende  niveaus  onderscheiden  worden. 


Tabel  1.1:  Definities  van  de  verschillende  niveaus  die  bij  doelopsporing  te  onderschei¬ 
den  zijn. 


detectie 

het  ontdekken  van  de  aanwezigheid  van  een  persoon,  object  of 
verschijnsel 

herkenning 

het  bepalen  van  de  algemene  aard  van  gedetecteerde  personen, 
obiecten  of  verschijnselen  (klasse  of  type) 

identificatie 

het  bepalen  van  de  specifieke  aard  van  gedetecteerde  personen, 
objecten  of  verschijnselen  (het  individu  binnen  een  klasse  of 

- - - 

In  het  algemeen  staat  doelopsporing  voor  detectie  en  herkenning.  Identificatie  als 
bijvoorbeeld  “vriend  of  vijand”  gaat  nog  een  stap  verder  en  vraagt  vaak  zeer  speci- 
fieke  kenmerken  en/of  een  vorm  van  communicatie.  Binnen  dit  onderzoek  is 
gekeken  naar  de  detectie  en  herkenning  van  stationaire  militaire  voertuigen. 

Het  onderzoek  bestond  uit  drie  fasen.  In  de  eerste  fase  is  een  cursus  gevolgd  op  het 
gebied  van  doeldetectie  en  herkenning  [Nasr  1996]  en  is  een  uitgebreide  literatuur- 
studie  gedaan  naar  bestaande  doelopsporingsmethoden  [van  Halsema  1998].  In  de 
tweede  fase  zijn  de  meest  belovende  methoden  geimplementeerd  en  toegepast  op 
PHARUS  data.  De  derde  fase  omvatte  de  evaluatie  van  de  toegepaste  methoden,  de 
rapportage  en  een  voorstel  voor  een  vervolgonderzoek.  Uit  de  eerste  fase  bleek  dat 
doelopsporing  als  volgt  weergegeven  kan  worden: 
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Figuur  1.1:  Doelopsporing  schematisch  weergegeven. 

De  laatste  drie  stappen  vormen  het  eigenlijke  doelopsporings  gedeelte.  In  de  pre¬ 
processing  worden  de  beelden  zoveel  mogelijk  ontdaan  van  mis.  Voor  SAR  is  dit 
een  belangrijke  stap  vanwege  de  sensor  inherente  speckle-mis.  In  de  detectie  wordt 
onderscheid  gemaakt  tussen  pixels  die  tot  een  doel  behoren  en  pixels  die  tot  de 
achtergrond  behoren.  Hierbij  zullen  ook  objecten  gedetecteerd  worden  die  niet 
relevant  zijn,  maar  dit  onderscheid  kan  worden  gemaakt  bij  de  herkenning.  Her- 
kenning  is  een  combinatie  van  feature  extractie  en  classificatie  op  basis  van  die 
features.  De  SAR  processing  en  calibratie  horen  niet  bij  doelopsporing  maar  zijn 
wel  weergegeven  omdat  zij  SAR  specifiek  zijn. 

Bij  de  indeling  van  doelopsporing  ligt  ook  de  scheiding  tussen  dit  onderzoek  en 
het  vervolgonderzoek.  Het  onderzoek  waar  hier  verslag  van  wordt  gedaan  omvat 
pre-processing  en  detectie.  Doelherkenning  wordt  kort  besproken  en  uitgewerkt. 
Het  vervolgonderzoek  gaat  hier  verder.  Door  een  vroege  afbakening  van  de  beide 
onderzoeken  (in  de  eerste  fase)  en  een  vertraging  van  het  hier  beschreven  onder¬ 
zoek  is  het  vervolgonderzoek  reeds  uitgevoerd  [de  Vries  1998].  Tussen  het  imple- 
menteren  van  de  meest  belovende  methoden  en  de  toepassing  op  PHARUS  data 
(tweede  fase),  zat  enige  tijd  omdat  de  calibratie  problemen  gaf  die  eerst  opgelost 
moesten  worden  [Dekker  1997-1,  1997-2]. 

Dit  rapport  volgt  ongeveer  de  lijn  van  figuur  1.1.  Het  eerst  volgende  hoofdstuk 
behandelt  het  polarimetrische  aspect  van  de  sensor.  In  hoofdstuk  3  wordt  de 
PHARUS  data  besproken.  Aan  de  orde  komen  de  (SAR)  processing,  de  calibratie, 
de  kwaliteit  van  de  data  en  de  ground-tmth.  Hoofdstuk  4  gaat  in  op  de  meest 
belovende  doelopsporingsmethoden  en  de  resultaten  hiervan  op  de  PHARUS  data. 
De  indeling  is  die  van  pre-processing,  detectie  en  herkenning.  In  hoofdstuk  5 
worden  conclusies  en  aanbevelingen  gedaan. 
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2.  Radar  polarimetrie 


Elektromagnetische  golven  hebben  een  orientatie  die  loodrecht  op  de  voortplan- 
tingsrichting  staat.  De  orientatie  richting  wordt  ook  wel  de  polarizatie  genoemd. 
Een  radar  die  zowel  horizontaal  als  verticaal  gepolarizeerde  golven  uitzendt  en 
ontvangt  wordt  een  polarimetrische  radar  genoemd.  Een  dergelijke  radar  kan 
maximaal  vier  kanalen  ontvangen:  HH,  HV,  VH  en  VV  (de  eerste  letter  staat  voor 
de  ontvangst  polarizatie,  de  tweede  letter  voor  de  zend  polarizatie).  De  samples  of 
pixels  kunnen  worden  beschreven  door  de  complexe  scattering  matrix: 


’vv  J 


A.,.,e 


JVi,. 


(1) 


A  is  de  gemeten  amplitude  en  (p  de  fase.  Voor  fase  gecalibreerde  data  geldt  dat 
Shy  =  Syh  [Zebker  1991].  De  samples  van  deze  data  kunnen  dus  ook  worden  be¬ 
schreven  door  de  complexe  vector: 


X  = 


‘-'/IV 

\SyyJ 


(2) 


Vermenigvuldigen  we  X  met  zijn  Hermitisch  getransponeerde,  dan  ontstaat  de 
covariantie  matrix: 


XX"" 


ShhShh 

ShyShh  ShySf,^ 
^SyyShh  ‘^VV‘^/1V 


ShhSyy 


ShySw 

SyySyy  ^ 


(3) 


H  staat  voor  Hermitisch  getransponeerd  en  *  voor  complex  geconjugeerd.  De 
Stokes  matrix  is  een  andere  notatie  van  de  covariantie  matrix.  Voor  de  meeste 
gedistribueerde  doelen  kan  (3)  ook  worden  geschreven  als  [Novak  1990]: 


^  1 
0 

y 


0  pVt' 
£  0 


(4) 


Hier  staat  Ohh  ^et  vermogen  van  HH,  e  voor  de  verhouding  tussen  het  ver- 
mogen  van  HH  en  HV,  yvoor  de  verhouding  tussen  het  vermogen  van  HH  en  W 
en  p  voor  de  correlatie  tussen  HH  en  VV.  De  resterende  producten  zijn  0  omdat  er 
doorgaans  geen  correlatie  bestaat  tussen  HH  en  HV  en  tussen  HV  en  VV.  Een 
voorwaarde  voor  deze  representatie  is  wel  dat  de  data  volledig  gecalibreerd  is 
[Dekker  1998].  Op  basis  van  de  parameters  van  (4)  is  onderscheid  te  maken  tussen 
doelen  van  verschillende  aard  of  type  (herkenning).  Hoe  goed  dit  gaat  is  een  deel 
van  de  gestelde  onderzoeksvraag  die  wordt  beantwoord  in  hoofdstuk  4. 
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3.  PHARUS  data  kwaliteit  en  ground-truth 


Voor  het  testen  van  de  doelopsporingsmethoden  die  in  het  volgende  hoofdstuk 
worden  behandeld  is  gekozen  voor  twee  PHARUS  opnamen  van  Swynnerton 
(UK).  Beiden  zijn  gemaakt  tijdens  de  NAVO  Panel  3  RSG-20  MIMEX  proefne- 
mingen  (30  September  t/m  25  oktober  1996).  Er  is  voor  deze  opnamen  gekozen 
omdat  zij  een  aantal  bekende  militaire  voertuigen  bevatten.  Daarbij  is  de  ground- 
truth  tijdens  de  proefnemingen  nauwkeurig  vastgelegd  [Bartlett  1997].  De  ruwe 
PHARUS  data  is  verwerkt  met  de  GSP  (Generieke  SAR  Processor)  [Otten  1994], 
een  ontwikkeling  van  TNO-FEL  en  ICT  Automatisering.  Vervolgens  zijn  de  op¬ 
namen  gecalibreerd  met  zelf  ontwikkelde  software  [Dekker  1998],  gebaseerd  op  de 
POLCAL  software  van  het  Jet  Propulsion  Laboratory  (USA)  [van  Zyl  1990]. 
Figuur  3.1  laat  de  opnamen  zien  van  het  gebied  waarin  de  voertuigen  zich  bevin- 
den  (ca.  750  m  x  750  m).  label  3.1  geeft  de  parameters  van  beide  opnamen  weer, 
na  de  verschillende  SAR  specifieke  operaties  (zie  figuur  1.1). 

Geometrie  en  signaal-ruis-vervouding 

Uit  figuur  3.1  en  tabel  3.1  blijkt  dat  de  opnamen  van  verschillende  kanten  zijn 
belicht.  Dit  verklaart  voor  het  grootste  gedeelte  de  verschillen.  De  opnamen  zijn 
met  een  tussentijd  van  circa  10  minuten  gemaakt.  Vooral  tijdens  opname  270X 
moest  de  antenne  flink  worden  bijgesuurd  (12°)  vanwege  de  heading  van  het 
vliegtuig  ten  opzichte  van  de  vliegrichting.  PHARUS  keek  bij  het  gebied  waarin  de 
doelen  zich  bevinden  (“area  of  interest”,  afgekort  a.o.i.)  zeer  scherend  in  (circa 
15°),  wat  de  lange  schaduwen  veroorzaakt.  Bij  deze  grazing-angle  is  de  minimale 
te  verwachten  gamma  -20  dB  (HV)  [Buiten  1990],  waaruit  valt  af  te  leiden  dat  de 
minimale  signaal-mis-verhouding  rond  de  1 1  dB  ligt  (zie  tabel  3.1).  Omdat  de 
gamma  van  HH  en  VV  gemiddeld  circa  7  dB  boven  die  van  HV  liggen,  volgt  dat 
de  minimale  signaal-ruis-verhouding  van  deze  kanalen  rond  de  18  dB  moet  liggen. 
Er  is  getracht  dit  na  te  gaan  voor  het  HH  kanaal,  door  de  gemiddelde  vermogens 
van  enkele  gedistribueerde  doelen,  te  vergelijken  met  die  binnen  enkele  schadu¬ 
wen,  zie  ook  figuur  3.1.  Schaduw  bevat  alleen  thermische  ruis  en  geen  speckle 
(multiplicatieve  ruis  die  het  gevolg  is  van  het  coherente  karakter  van  SAR  data, 
normaal  gesproken  sterker  dan  thermische  mis).  Hiemit  bleek  dat  deze  verhouding 
10  dB  slechter  was  (8  dB).  Blijkbaar  bevatten  de  schaduwen  niet  alleen  thermische 
mis.  Uit  een  nadere  analyse  van  de  data  bleek  ten  eerste  dat  de  mwe  data  was 
“vastgelopen”.  De  versterkers  van  PHARUS  worden  zo  ingesteld  dat  de  mis  van 
de  ontvangers  net  boven  de  kwantisatie-mis  van  de  AD-converters  (ADC)  uitkomt 
[Vermeulen  1996].  Voor  de  zekerheid  was  hier  10  dB  extra  bij  opgeteld.  Bij 
gebmik  van  de  “most-significant”  8  bits  van  de  data  gaat  daardoor  niets  verloren. 
De  Gaussisch-vormige  histogrammen  van  de  I  en  Q  kanalen  waren  echter  zo  breed 
(standaard  deviaties  van  rond  de  90)  dat  zij  abmpt  worden  afgesneden  bij  ±127. 
Wat  de  gevolgen  hiervan  zijn  is  moeilijk  na  te  gaan,  maar  dat  het  zorgt  voor  een 
slechter  beeld  is  te  zien  in  gebieden  waar  de  reflecties  hoger  zijn  (gebieden  met 
bijvoorbeeld  gebouwen).  De  achtergronden  van  deze  reflecties  worden  daarbij 
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270X 


vliegrichting  - 1 

I  kijkrichting 

Figuur  3.1:  De  PHARUS  opnamen  van  Swynnerton  (UK)  waarmee  de  doelopsporings- 
methoden  warden  getest  (R=HH,  G=HV,  B=W).  In  het  midden  bevindt  zich 
het  gebied  waar  de  voertuigen  (de  doelen)  zijn  opgesteld. 
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Tabel  3. 1:  De  (kwaliteit)  parameters  van  de  PHARUS  opnamen  van  Swynnerton  ( UK) 

tijdens  de  verschillende  SAR  specifieke  operaties. 


common 

180X 

270X 

date 

d/m/y 

24/10/1997 

time 

GMT 

12:56:36 

12:45:54 

track 

O 

O 

00 

270° 

beam  angle 

-4.5° 

O 

p 

altitude 

3096  m 

3052  m 

projection 

slant  range 

slant  range 

a.o.i. 

12190  m 

11490  m 

grazing  angle 

a.o.i. 

14.7° 

15.4° 

signal-to-noise 

Y=-20  dB 

10.4  dB 

11.6  dB 

accuracy 

ADC 

8  bits 

deviation 

HH 

97.4 

88.0 

after  processing 

pixel  spacing 

range 

1.5  m 

azimuth 

0.31518  m 

0.35770  m 

resolution 

range 

5.0  m 

azimuth 

1.2  m 

ENL 

HH 

1 

a.o.i. 

175° 

a.o.i. 

135° 

cross-talk  level 

<  -20  dB 

gain-imbalance 

-5.3  dB 

-3.5  dB 

after  calibration 

pixel  spacing 

range 

1.5  m 

azimuth 

1.5759  m 

1.4308  m 

resolution 

range 

5.0  m 

azimuth 

3.1518  m 

3.5770  m 

ENL 

HH 

2.8 

(Pu,  -  <Pvv 

<  10° 

cross-talk  level 

<  -30  dB 

gain-imbalance 

<  1.0  dB 
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sterk  onderdrukt.  In  gebieden  waar  dit  minder  goed  is  te  zien  kan  het  echter  nog 
steeds  zo  zijn  dat  de  gevolgen  van  het  vastlopen  van  de  ruwe  data,  de  oorspronke- 
lijke  signaal-ruis-verhouding  overtreft.  Ten  tweede  was  er  tijdens  de  opnamen 
geen  tapering  toegepast  waardoor  de  zijlussen  van  de  antenne  niet  zijn  onderdrukt. 
Dit  kan  tot  gevolg  hebben  dat  de  zijlussen  van  een  doel  doorspreken  in  een  ander 
doel.  Tot  slot  spreken  ook  de  randen  door  en  zijn  de  schaduwen  mogelijk  te  kort 
(tot  50  m  slant  range),  om  een  goede  bepaling  te  doen  van  de  signaal-ruis- 
verhouding 

SAR  processing 

De  ruwe  PHARUS  data  ondergaat  twee  correlatie  slagen  voordat  het  een  scherp 
beeld  oplevert.  Er  is  gekozen  voor  een  maximale  resolutie  met  weging:  1.0  m  in 
azimuth  (along  track)  en  3.4  m  in  range  (cross-track)  richting.  De  oorspronkelijke 
pixel  afmeting  in  range  van  de  complexe  ruwe  PHARUS  data  bedroeg  3.0  m.  Na 
de  correlatie  en  calibratie  is  de  data  echter  niet  meer  complex  maar  gedetecteerd 
zodat  dit  niet  meer  voldoende  is  ten  opzichte  van  de  resolutie.  Daarom  is  de  pixel 
afmeting  van  de  complexe  ruwe.  data  (en  de  chirp)  gehalveerd  tot  1.5  m  door 
interpolatie.  Gekozen  is  voor  de  nearest-neighbour  methode  omdat  deze  de  statisti- 
sche  eigenschappen  (gemiddelde,  standaard  deviatie,  autocorrelatie  funtie,  etc.)  het 
beste  behoudt  [Quegan  1990].  Na  metingen  bleek  de  range  resolutie  geen  3.4  m 
maar  5.0  m.  De  oorzaak  hiervan  bleek  het  gedeelte  van  de  software  dat  zorgt  voor 
de  range-migratie  correctie  tijdens  de  azimuth  correlatie.  V oor  de  azimuth  correla- 
tie  bleek  de  range  resolutie  tussen  3.0  m  en  4.0  m  te  liggen.  De  pixel  afmeting  in 
azimuth  bleek  ruim  voldoende  voor  een  azimuth  resolutie  die  1.2  m  bleek  te  be- 
dragen. 

Calibratie 

De  calibratie  van  de  beelden  houdt  het  volgende  in: 

•  fase  calibratie 

•  receiver-gain  correctie 

•  symmetrisatie 

•  conversie  naar  covariantie  matrix 

•  multi-looking  (optioneel) 

•  cross-talk  calibratie 

•  gain-imbalance  calibratie 

•  absolute  calibratie 

Voor  details  wordt  verwezen  naar  [Dekker  1998].  De  eerste  stap  trekt  de  fases  van 
de  vier  kanalen  in  de  scattering  matrix  gelijk,  zie  vergelijking  (1).  De  receiver-gain 
correctie  compenseert  voor  de  instelling  van  de  versterkers  die  vooraf  gaan  aan  de 
AD-converters  van  PHARUS.  In  de  symmetrisatie  worden  de  kanalen  HV  en  VH 
gemiddeld,  ter  verbetering  van  de  signaal-ruis-verhouding  in  het  nieuwe  HV 
kanaal.  Multi-looking  is  optioneel  en  houdt  in  het  middelen  van  een  aantal  samples 
in  azimuth  zodat  de  speckle  afneemt.  Door  multi-looking  neemt  ook  de  resolutie  af 
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en  kan  het  data  volume  worden  verkleind.  In  de  volgende  paragraaf  wordt  dieper 
ingegaan  op  multi-looking.  De  cross-talk  calibratie  zuivert  de  verzonden  en  ont- 
vangen  polarizaties  van  andere  orientaties.  Uitgezonden  horizontaal  gepolariseerde 
golven  bijvoorbeeld,  blijken  ook  een  kleine  verticale  component  te  hebben.  De 
gain- imbalance  calibratie  corrigeert  voor  gain  verschillen  tussen  het  horizontale  en 
verticale  pad  binnen  de  radar,  meestal  veroorzaakt  door  verschillen  in  antenne-gain 
voor  H  en  V.  Bij  PHARUS  is  de  gain-imbalance  erg  groot  en  afhankelijk  van  de 
stand  van  de  antenne-bundel  (beam  angle).  Dit  komt  omdat  een  andere  stand  van 
de  antenne-bundel  resulteert  in  een  ander  antennepatroon  [Vermeulen  1996]  en 
dus  een  andere  gain.  De  gain-imbalance  van  de  hier  besproken  opnamen  bedroeg  - 
5.3  dB  en  -3.5  dB  voor  calibratie,  corresponderend  met  beam  angles  van  -4.5°  en 
12°.  In  de  absolute  calibratie  wordt  het  ontvangen  vermogen  omgerekend  naar 
gamma  waarden  (van  den  Broek  1993].  De  absolute  calibratie  is  niet  uitgevoerd 
omdat  het  absolute  nivo  niet  van  belang  is  binnen  de  covariantie  matrix  gegeven 
door  vergelijking  (4). 


Eauivalent-number-of-looks 

Ten  behoeve  van  de  presentatie  van  de  data  is  gekozen  voor  een  vierkante  resolutie 
cel  van  ongeveer  3.4  m  x  3.4  m  (dat  de  range  resolutie  5.0  m  bedroeg  bleek  later). 
Hierdoor  was  het  mogelijk  om  iedere  10  samples  in  azimuth  te  middelen  en  een 
reductie  factor  van  4  a  5  toe  te  passen.  Dit  wordt  ook  wel  multi-looking  genoemd. 
Het  equivalent-number-of-looks  (ENL)  is  een  maat  voor  de  speckle  reductie  die 
hierbij  optreedt,  en  wordt  gedefmieerd  als  [Moreira  1991]: 


^  is  het  gemiddelde  en  crde  standaard  deviatie.  Hierbij  geldt  dat  het  gebied  waar- 
binnen  deze  parameters  worden  bepaald,  homogeen  moet  zijn  en  geen  textuur  mag 
bevatten.  Of  dit  het  geval  is  zal  proefondervindelijk  moeten  worden  vastgesteld, 
door  het  ENL  te  bepalen  van  meerdere  ogenschijnlijk  homogene  gebieden  (wat  op 
het  oog  homogeen  is  hoeft  nog  niet  zo  zijn).  De  meest  waarschijnlijke  ENL  zal 
zich  bevinden  bij  de  maximale  waarden.  De  waarde  2.8  in  tabel  3.1  (na  calibratie) 
is  met  deze  methode  bepaald.  Het  ENL  kan  ook  worden  berekend  met  behulp  van 
de  vergelijking  [Moreira  1991]: 


ENL  = 


1  +  2IVp;.u 

;=i  ^ 


(6) 


L  is  het  aantal  samples  dat  wordt  gemiddeld  (hier  10)  en  is  de  correlatie 
coefficient  tussen  sample  i+1  en  sample  1.  De  correlatie  coefficient  kon  worden 


afgeleid  met  behulp  van  [Ulaby  1986]: 


Pma  =  sinc- 


.  J 


(7) 
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Hier  is  dx  de  pixel  afstand  in  azimuth  en  de  azimuth  resolutie.  Bij  de  hier 
toegepaste  multi-look  methode  is  alleen  de  correlatie  in  azimuth  richting  van 
belang.  Door  de  gegevens  van  tabel  3.1  in  te  vullen  in  vergelijking  (6)  en  (7) 
volgen  een  ENL  van  4.3  voor  180X  en  4.8  voor  270X.  De  waarde  bepaald  met  de 
methode  van  vergelijking  (5)  is  echter  significant  lager;  waarschijnlijk  veroorzaakt 
door  de  slechtere  signaal-ruis-verhouding  ten  gevolge  van  het  vastlopen  van  de 
data  en  het  niet  toepassen  van  antenne  tapering.  Het  is  niet  waarschijnlijk  dat  alle 
gebieden  waarbinnen  het  ENL  werd  bepaald,  minder  homogeen  waren  dan  ver- 
wacht.  Door  bij  het  middelen  van  samples  gebruik  te  maken  van  weging,  is  het 
mogelijk  om  (met  behoud  van  resolutie)  het  ENL  te  verhogen  met  ongeveer  een 
factor  1.5  [Li  1983].  Dit  is  niet  verder  onderzocht  omdat  de  tijd  hiervoor  ontbrak. 

Ground-truth 

Binnen  dit  onderzoek  is  gekeken  naar  de  waameembaarheid  van  stationaire  mili- 
taire  voertuigen  met  PHARUS.  Hierbij  is  data  gebruikt  die  is  opgenomen  tijdens 
de  in  de  inleiding  met  naam  genoemde  proefnemingen  bij  Swynnerton  (UK).  Een 
van  de  redenen  waarom  gekozen  is  voor  deze  opnamen  is  de  aanwezigheid  van 
ground-truth  data  [Bartlett  1997].  Tijdens  deze  proefnemingen  zijn  er  in  een  ter- 
rein  van  ca.  250  m  x  250  m  vijf  militaire  voertuigen  neergezet.  Het  terrein  betrof 
bos  afgewisseld  door  veel  grote  en  kleine  open  plekken  hoofdzakelijk  begroeid 
met  gras.  Op  dit  terrein  bevond  zich  ook  een  baan  voor  het  werpen  van  handgra- 
naten  (grenade-range).  Bijlage  A  toont  een  collectie  foto’s  en  een  kaart  van  het 
terrein.  De  collectie  omvat  luchtfoto’s,  panorama’s  en  foto’s  van  de  gestationeerde 
voertuigen.  De  locaties  van  de  voertuigen  zijn  aangegeven  op  de  kaart.  Op  basis 
van  deze  informatie  is  in  uitsneden  van  de  twee  PHARUS  opnamen  van  figuur  3.1 
aangegeven  wat  de  voertuigen  zijn,  zie  figuur  3.2.  Hierbij  valt  op  dat  dit  onder- 
scheid  met  het  oog  niet  gemakkelijk  is  te  maken.  In  het  volgende  hoofdstuk  wordt 
bekeken  of  deze  doelen  gedetecteerd  kunnen  worden  met  behulp  van  algoritmen  en 
of  zij  te  onderscheiden  zijn  van  andere  objecten  die  geen  doel  zijn  maar  wel  gede¬ 
tecteerd  worden  (false-alarms),  anders  gezegd:  of  zij  herkend  worden  als  doel. 
Door  de  resolutie  van  de  PHARUS  opnamen  is  het  aantal  pixels  dat  een  doel 
beslaat  klein.  Dit  sluit  echter  niet  bij  voorbaat  uit  dat  PHARUS  ongeschikt  is  voor 
het  opsporen  van  doelen.  Doeldetectie  en  herkenning  zijn  vaak  een  kwestie  van  het 
gebruiken  van  meerdere  sensoren  (sensor-fusie),  waarbij  sensoren  elkaar  onder- 
steunen  [Nasr  1996].  Een  voorbeeld  hiervan  is  dat  een  sensor  een  doel  detecteert 
en  aangeeft  aan  andere  sensoren  ter  herkenning  (cueing). 
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180X 


270X 


1  =  GECKO 

2  =  TEAR 

3  =  ASTROSS 

4  =  MTLB 

5  =  ZSU  23-4 

Figuur  3.2:  Uitvergroting  (3x)  van  de  PHARUS  opnamen  vanfiguur  3. 1  waarin  de 
voertuigen  zijn  aangegeven. 
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4.  Doelopsporing  theorie  en  resultaten 

In  dit  hoofdstuk  worden  de  meest  belovende  doelopsporingsmethoden  besproken 
en  toegepast  op  de  besproken  PHARUS  opnamen.  De  onderverdeling  is  die  van 
figuur  l.l:  pre-processing  methoden,  detectoren  en  herkenning. 


4.1  Pre-processing:  speckle  jRItering 

Detectie  in  SAR  beelden  heeft  baat  bij  de  reductie  van  speckle,  omdat  de  stan- 
daarddeviatie  van  de  achtergrond  van  een  doel  hierdoor  afneemt  (zie  ook  paragraaf 
4.2).  Binnen  dit  onderzoek  is  een  literatuur  studie  gedaan  naar  polarimetrische  en 
spatieel-polarimetrische  filters.  Het  verslag  hiervan  is  te  vinden  in  [van  Halsema 
1998].  De  meest  belovende  methoden  zijn  geimplementeerd  en  worden  in  de 
volgende  alinea’s  beschreven.  Belangrijk  bij  filteren  is  dat  de  resolutie  gelijk  blijft 
zodat  de  reflecties  van  de  doelen  niet  aangetast  worden.  Daarom  zijn  bijna  alle 
spatiele  methoden  afgevallen.  De  hier  beschreven  filters  zijn  dan  ook  allemaal 
polarimetrische  filters,  op  het  wavelet-  shrinkage  filter  na.  Dit  filter  bleek  zulke 
goede  resultaten  op  te  leveren  [Guo  1994]  en  was  dermate  onbekend,  dat  het  is 
meegenomen  in  dit  onderzoek. 

Standard-deviation-to-mean  ratio 

De  filters  zullen  worden  vergeleken  aan  de  hand  van  de  standard-deviation-to- 
mean  ratio  die  veel  gebruikt  wordt  in  de  bestudeerde  literatuur.  De  standard-devia¬ 
tion-to-mean  wordt  gedefinieerd  als  [Novak  1990]: 

s/m  =  — 

Vergelijken  we  deze  definitie  met  die  van  vergelijking  (5)  dan  volgt  dat  de  stan¬ 
dard-deviation-to-mean  gelijk  is  aan  I/VENL  .  Net  als  bij  vergelijking  (5)  geldt 
ook  hier  dat  het  gebied  waarbinnen  de  standaard  deviatie  en  het  gemiddelde  wor¬ 
den  bepaald,  homogeen  moet  zijn  en  geen  textuur  mag  bevatten.  De  standard- 
deviation-to-mean  ratio,  s/m  ratio  of  A(s/m)  is  het  verschil  tussen  de  in-  en  output 
standard-deviation-to-mean  van  een  filter.  Beiden  (s/m  en  ratio)  worden  meestal 
uitgedrukt  in  dB’s  (201og(s/m)). 

Dual  Channel 

Polarimetrische  data  is  normaal  gesproken  driekanaals.  Er  zijn  ook  systemen  te 
bedenken  die  twee  kanalen  bezitten.  In  [Nasr  1990]  is  onderzocht  wat  dit  voor 
effect  heeft  op  het  detecteren  van  gedistribueerde  doelen.  Onderzocht  is  het  ge¬ 
middelde: 

|////|'  +  |VTp 


2 


(9) 
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Voor  de  standard-deviation-to-mean  ratio  geldt: 

(.0) 

'  1  +  y 

Voor  de  parameters  yen  p  wordt  verwezen  naar  hoofdstuk  2.  Het  is  gemakkelijk  in 
te  zien  dat  wanneer  beide  kanalen  gelijk  zijn  (  y  =  1 )  en  niet  gecorreleerd  ( p  =  0 ), 
de  s/m  ratio  maximaal  -3  dB  kan  worden.  In  de  praktijk  komen  dergelijke  achter- 
gronden  zelden  voor  en  is  de  s/m  ratio  meestal  slechter  dan  -3  dB. 


Span 

Het  gemiddelde  vermogen  van  driekanaals  data,  ook  wel  de  span  of  total  power 
genoemd,  wordt  doorgaans  als  volgt  berekend  [Nasr  1990,  Novak  1990]: 

_  !////!-  +  2\HV[  +  \W[  ( 

4 


Hierbij  wordt  het  vermogen  van  HV  twee  keer  meegenomen  omdat  dit  gemiddeld 
lager  is  dan  dat  van  HH  en  VV.  Voor  de  s/m  ratio  geldt: 


A(s/m)  = 


.^l  +  4£^+y^+2|ppy 

l  +  2e  +  y 


(12) 


De  maximale  s/m  ratio  voor  het  span  filter  volgt  voor  £  =  0.5 ,  y  =  1  en  p  =  0 ,  en 
bedraagt  -4.77  dB.  Hierbij  mag  er  geen  correlatie  zijn  tussen  HV  en  de  overige 
twee  kanalen.  In  de  praktijk  is  de  s/m  ratio  bijna  altijd  slechter  dan  -4.77  dB. 


OIS 

Het  Optimal  Intensity  Sum  (OIS)  filter  is  vergelijkbaar  met  de  vorige  twee  filters. 
Het  verschil  is  echter  dat  het  een  adaptief  filter  is  en  een  optimaal  gemiddelde  of 
som  van  de  vermogens  berekent  op  basis  van  de  covariantie  matrix  van  de  gemid¬ 
delde  achtergrond.  De  filtervergelijking  is  als  volgt  [Novak  1990]: 


\HHf  + 
y  = - 


+  -1  wp 

£ _ y 

3+ipr 


(13) 


Dit  filter  kent  twee  varianten:  (a)  waarbij  in  een  moving  window  rond  ieder  pixel 
de  gemiddelde  covariantie  matrix  van  de  achtergrond  wordt  berekend  en  (b)  een 
snelle  variant  waarbij  een  gemiddelde  covariantie  matrix  wordt  gebruikt  voor  het 
hele  beeld.  Voor  de  s/m  ratio  van  de  (a)  variant  geldt: 


V3+4ipr+ipr 

yj9  +  6\pf  +\p\' 


(14) 


Voor  de  (b)  variant  geldt  deze  vergelijking  alleen  als  de  gemiddelde  covariantie 
matrix  rondom  een  pixel  gelijk  is  aan  de  globale  covariantie  matrix.  De  voorkeur 
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gaat  uit  naar  de  (b)  variant  vanwege  de  snelheid  en  het  feit  dat  de  te  detecteren 
doelen  in  een  beeld  vaak  eenzelfde  achtergrond  hebben.  Is  er  geen  correlatie  tussen 
HH  en  VV  ( p  =  0 )  dan  is  de  s/m  ratio  -4.77  dB.  Zijn  zij  volledig  gecorreleerd 
( p  =  1 )  dan  is  de  s/m  ratio  -3  dB. 


PWF 

Het  Polarimetric  Whitening  Filter  [Novak  1990,  1993]  is  het  optimale  polarimetri- 
sche  filter.  De  essentie  van  dit  filter  is  dat  het  de  oorspronkelijke  polarimetrische 
data  met  de  lineaire  basis  (HH,  HV,  VV)  transformeert  naar  data  met  een  basis 
waarin  de  drie  kanalen  niet  gecorreleerd  zijn.  Het  PWF  beeld  is  het  gemiddelde 
van  deze  kanalen.  De  filtervergelijking  voor  multi-look  data  wordt  gegeven  door 
[Liu  1995]: 

»;(X^  (15- 

3 


Hier  is  XX"  de  sample  covariantie  matrix,  zie  (3),  Xc  de  gemiddelde  covariantie 
matrix  van  de  omgeving  rond  XX''  en  tr{ )  de  trace  van  een  matrix.  Het  PWF  filter 
kent  dezelfde  twee  varianten  als  het  OIS  filter.  Voor  de  s/m  ratio  van  de  (a)  variant 
geldt: 


A(s/m)  =  ^  (-4.77  dB) 


(16) 


Om  dezelfde  redenen  als  bij  het  OIS  filter  gaat  de  voorkeur  hier  ook  uit  naar  de  (b) 
variant.  Het  verschil  in  performance  tussen  beide  varianten  bleek  minimaal  [Novak 
1993].  Vergelijking  (16)  geldt  dus  ook  voor  de  (b)  variant.  Voor  de  covariantie 
matrix  Zc  kan  het  beste  worden  gekozen  voor  die  van  de  achtergrond  van  de  te 
detecteren  doelen. 


DIWF 

Polarimetrische  data  kan  ook  getransformeerd  worden  naar  data  met  de  lineaire 
basis  (HH-VV,  HV,  HH+VV).  De  eerste  twee  kanalen  van  deze  basis  behoren  tot 
de  dihedral-comer  (even)  reflecties,  het  laatste  kanaal  tot  de  vlakke  plaat  en  trihe¬ 
dral-comer  (oneven)  reflecties  [Freeman  1992].  Hierbij  slaat  even  en  oneven  op 
het  aantal  keren  dat  de  radar  golven  worden  weerkaatst.  Dihedral-  en  trihedral- 
comer  reflecties  zijn  kenmerkend  voor  man-made  objecten  als  militaire  voertui- 
gen.  Analoog  aan  het  PWF  filter  wordt  de  optimale  combinatie  van  de  kanalen 
(HH-W,  HV)  het  Dihedral  Whitening  Filter  (DIWF)  genoemd.  Het  filter  kan  als 
volgt  worden  geschreven  [Larson  1995]: 

_ 


l  +  y-2VyRe(p) 


De  s/m  ratio  van  dit  filter  is  niet  bekend. 
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TRIWF 

Analoog  aan  het  DIWF  is  ook  het  Trihedral  Whitening  Filter  (TRIWF)  uitgewerkt: 

3;  = - !— - \HH  +  Wf  (18) 

l  +  Y  +  2.ylyRs{p) 

De  s/m  ratio  van  dit  filter  is  net  als  bij  het  DIWF  niet  bekend. 

Wavelet  Shrinkage  filter 

Het  enige  niet  polarimetrische  filter  dat  is  geimplementeerd  is  het  Wavelet  Shrin¬ 
kage  of  VisuShrink  filter  [Donoho  1994,  1995].  Het  filter  is  uitgewerkt  ten  behoe- 
ve  van  doel  detectie  en  herkenning  [Guo  1994].  Het  Wavelet  Shrinkage  filter  komt 
in  essentie  neer  op  “soft-thresholding”  in  het  wavelet  domein,  zie  figuur  4.1. 


Figuur  4. 1:  Schematische  weergave  van  het  Wavelet  Shrinkage  filter . 

Hierbij  wordt  de  wavelet  transformatie  uitgevoerd  op  logaritmisch  geschaalde 
data.  Dit  heeft  het  voordeel  dat  de  speckle  die  van  nature  multiplicatief  is,  additief 
wordt.  Nadeel  hiervan  is  dat  het  gemiddelde  na  de  inverse  logaritme  niet  meer  te 
reconstrueren  is  [Hoekman  1991].  Voor  detectie  is  dit  geen  probleem.  Andere 
operaties  als  classificatie  zouden  hierdoor  echter  wel  verkeerd  kunnen  gaan. 

Als  wavelet  transformatie  werd  gekozen  voor  de  Daubechie’s  wavelet  met  vier 
coefficienten  (D4)  [Press  1992].  Een  wavelet  transformatie  is  de  inverse  Fourier 
transformatie  van  een  Fourier  (frequentie)  domein  dat  is  ingedeeld  in  verschillende 
logaritmisch  ingedeelde  frequentie  intervallen  [Chui  1992].  Zij  bewerkstelligt  een 
decompositie  van  de  verschillende  resolutie  niveaus  die  in  een  beeld  aanwezig 
zijn,  zie  figuur  4.2.  Het  hoogste  niveau  (n-1)  bevat  voomamelijk  ruis  en  het  laagste 
niveau  (0)  is  de  DC  component;  het  gemiddelde.  Alleen  de  hoogste  niveaus  (op 
zijn  minst  vijf)  ondergaan  soft-thresholding.  De  laagste  niveaus  (vanaf  0)  blijven 
gelijk  waardoor  het  gemiddelde  van  het  beeld  behouden  blijft.  Bij  soft- 
thresholding  worden  waarden  niet  begrensd,  zoals  bij  hard-thresholding  (clipping), 
maar  verlaagd  met  een  bepaalde  waarde: 

jsign(w)(lwl-/)  iflvvl>/ 

“I  0  iflwl</ 

Hier  is  t  de  soft-threshold  en  vv  een  sample  van  een  van  de  hoogste  niveaus  in  het 
wavelet  domein.  [Guo  1994]  gaf  aan  dat  de  threshold  t  het  beste  kon  liggen  tussen 
de  l.Scren  ScTomdat  de  doelen  anders  te  veel  werden  aangetast.  De  standaarddevi- 
atie  werd  bepaald  in  het  hoogste  niveau  (n-1).  De  s/m  ratio  van  dit  filter  laat  zich 
moeilijk  voorspellen.  In  de  literatuur  [Guo  1994]  worden  met  de  beschreven  me- 
thode  s/m  ratios  behaald  van  circa  0.5  (-6  dB)  bovenop  PWF  gefilterde  data. 
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Figuur  4.2: 


n-2  n-1 


n-2 


n-1 


Het  wavelet  domein  van  een  tweedimensionale  wavelet  transformatie. 


Resultaten 

Alle  besproken  filters  behalve  het  dual  channel  en  span  filter  zijn  adaptief  Het 
OIS  en  de  whitening  filters  (PWF,  DIWF  en  TRIWF)  hebben  de  covariantie  matrix 
van  de  achtergrond  nodig  als  input.  In  beide  PHARUS  opnamen  (zie  hoofdstuk  3) 
is  nagegaan  wat  de  gemiddelde  (20  x  20  pixels)  achtergrond  matrix  was,  zie  tabel 
4.1,  kolom  Ec-  Uit  de  foto’s  van  bijlage  A  bleek  de  achtergrond  gras  te  zijn.  De  s/m 
ratio’s  in  dezelfde  kolom  zijn  berekend  met  de  formules  (10),  (12),  (14)  en  (16)  en 
de  covariantie  matrix.  In  deze  tabel  zijn  ook  weergegeven  de  covariantie  matrices 
van  enkele  homogene  testgebieden  die  zijn  gebruikt  voor  het  bepalen  van  de  s/m 
ratio’s  van  de  filters.  Aan  de  covariantie  matrices  te  zien  zijn  de  velden  waar- 
schijnlijk  verschillend,  de  exacte  aard  is  echter  niet  bekend.  Het  Wavelet  Shrinka¬ 
ge  filter  is  ook  adaptief  en  bepaalt  zijn  resultaat  op  basis  van  de  standaard  deviatie 
van  de  speckle.  Het  is  gei'mplementeerd  als  beschreven.  De  lengte  van  de  D4 
wavelet  bedroeg  512  x  512  samples.  Soft-thresholding  werd  toegepast  op  de 
hoogste  vijf  niveaus  (4  t/m  8).  De  standaard  deviatie  werd  bepaald  op  het  hoogste 
niveau  (8).  Uitgewerkt  zijn  de  fitters  met  een  soft-threshold  van  2.5cr  (W2.5S)  en 
5a  (W5S). 


Bijlage  B  toont  het  HH,  HV  en  W  kanaal  van  de  opnamen  en  de  resultaten  van 
alle  hier  beschreven  filters.  Tabel  4.1  toont  de  corresponderende  s/m  ratio’s  ge- 


meten  in  de  testgebieden.  We  zien  dat  de  polarimetrische  filters  binnen  een  kleine 
marge  voldoen  aan  de  verwachtingen.  Verschillen  kunnen  worden  veroorzaakt 
doordat  er  toch  enige  textuur  aanwezig  is,  of  door  de  slechtere  signaal-ruis- 
verhouding  besproken  in  hoofdstuk  3.  De  verschillen  blijven  echter  klein.  Dit  geldt 
niet  voor  het  DIWF  en  het  TRIWF.  De  Wavelet  Shrinkage  filters  presteren  minder 
dan  die  van  [Guo  1994],  Vermoed  wordt  dat  de  prestaties  van  deze  filters  afhangen 
van  wat  er  aan  speckle  wordt  aangeboden,  ofwel,  dat  wanneer  er  minder  mis 
aanwezig  is,  er  minder  te  filteren  valt.  In  beide  gevallen  werd  PWF  data  aangebo¬ 
den,  echter  het  PWF  hier  gaf  betere  resultaten  dan  bij  [Guo  1994],  Vergelijken  we 
onze  resultaten  van  beide  wavelet  filters  met  die  van  het  PWF  (bijlage  B),  dan  zien 
we  dat  het  2.5a  filter  een  lichte  verbetering  geeft  en  dat  het  5ct  filter  het  beeld 
begint  aan  te  tasten.  De  grens  van  3(7  lijkt  daardoor  bevestigd. 

Uit  tabel  4.1  blijkt  dat  het  PWF  filter  het  beste  polarimetrische  filter  is.  Span  en 
OIS  zijn  hiervoor  een  goed  altematief,  mede  door  de  eenvoud  van  deze  filters.  Het 
Wavelet  Shrinkage  filter  geeft  nog  enige  verbetering  op  het  PWF.  Of  het  Wavelet 
Shrinkage  filter  doelen  aantast  zal  moeten  blijken  bij  de  detectie  in  paragraaf  4.2. 
Het  DIWF  presteert  goed,  vooral  bij  bomen.  Het  TRIWF  blijkt  helemaal  niet  te 
filteren.  De  combinatie  DIWF/TRIWF  kan  worden  gezien  als  een  decompositie 
van  dihedral  (even)  reflecties  en  vlakke  plaat  en  trihedral  (oneven)  reflecties. 
Figuur  4.3  toont  een  compositie  van  beide  beelden  in  complementaire  kleuren.  Een 
dergelijk  beeld  geeft  een  indicatie  of  reflecties  meer  even  zijn  dan  oneven  of 
omgekeerd.  Wat  opvalt  is  dat  vooral  in  de  PHARUS  opname  270X  er  enkele 
voertuigen  zijn  met  een  uitgesproken  aard  (even/oneven).  Tabel  4.2  laat  tot  slot 
zien  wat  de  equivalent-number-of-looks  zijn  van  de  oorspronkelijke  data  en  na 
toepassing  van  de  verschillende  filters.  Wat  opvalt  is  dat  de  equivalent-number-of- 
looks  van  de  kanalen  HH,  HV  en  VV  niet  gelijk  zijn.  Dit  kan  doorspreken  in  de 
resultaten  van  de  diverse  filters. 
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Tabel  4.1:  De  covariantie  matrices  van  de  achtergrond  Z,.  en  de  test  gebieden,  en  de 

s/m  ratio's  van  deze  gebieden  na  de  verschillende  filters. 


2c 

180X 

field  trees 

field 

270X 

trees 

£ 

0.19 

0.16 

0.32 

0.10 

0.27 

Y 

1.00 

1.64 

0.93 

0.92 

0.86 

\P\ 

0.45 

0.62 

0.36 

0.38 

0.29 

0° 

2.7° 

26.1° 

23.9° 

28.3° 

s/m  to  HH 

TRIWP 

-0.11 

1.38 

0.72 

1.13 

DUAL 

-2.21 

-1.86 

-1.39 

-2.10 

-1.46 

DIWF 

-2.77 

-4.41 

-2.22 

-4.01 

SPAN 

-3.47 

-2.84 

-2.78 

-2.86 

-3.18 

OIS 

-4.25 

-3.67 

-3.82 

-3.71 

-4.04 

PWF 

-4.77 

-A.n 

-4.82 

-4.23 

-4.61 

s/m  to  PWF 

W2.5S 

-1.74 

-1.16 

-1.02 

-1.51 

W5S 

-2.93 

-1.83 

-1.66 

-2.82 

Tabel  4.2:  D 

fi 

e  equivalent-number-of-looks  van  de  test  gebieden  na  de  verschillende 
Iters. 

180X 

field  trees 

270X 

field  trees 

HH 

2.74 

2.20 

3.01 

2.84 

HV 

1.72 

2.49 

2.50 

3.43 

VV 

3.16 

2.00 

2.04 

2.49 

TRIWF 

2.81 

1.60 

2.55 

2.19 

DUAL 

4.21 

3.02 

4.88 

3.97 

DIWF 

5.18 

6.08 

5.02 

7.15 

SPAN 

5.26 

4.17 

5.81 

5.90 

OIS 

6.37 

5.29 

7.07 

7.19 

PWF 

8.20 

6.66 

7.97 

8.20 

PWF+W2.5S 

PWF+W5S 


12.24 

16.11 


8.69 

10.15 


10.07 

11.67 


11.61 

15.69 
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180X 


270X 


1  =  GECKO 

2  =  TEAR 

3  =  ASTROSS 

4  =  MTLB 

5  =  ZSU  23-4 


Figuur  4.3:  Composities  van  de  DIWF  beelden  ( rood)  en  TRIWF  beelden  ( cyaan):  twee 
decomposities  van  dihedral  (even)  reflecties  en  trihedral  (oneven)  reflecties. 
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4.2  Detectoren 

Na  de  pre-processing  (speckle  filtering)  volgt  de  detectie  van  de  samples  (pixels) 
die  tot  een  doel  behoren.  Verschillende  detectoren  zijn  onderzocht. 


Klassieke  CFAR  detector 

Bij  veel  detectoren  wordt  de  waarde  van  een  sample  vergeleken  met  het  gemiddel- 
de  van  de  achtergrond  clutter.  Is  het  verschil  groter  dan  een  bepaalde  threshold  dan 
wordt  aangenomen  dat  deze  sample  tot  een  doel  behoort.  Bij  Constant  False  Alarm 
Rate  (CFAR)  detectoren  is  de  threshold  gesplitst  in  de  standaard  deviatie  van  de 
achtergrond  en  een  constante  [Goldstein  1973]: 


^  ^CFAR 


(20) 


Hier  is  x,  de  sample  die  getest  wordt,  ^  het  gemiddelde  van  de  achtergrond,  Oc  de 
standaard  deviatie  van  de  achtergrond  en  Kcfar  de  CFAR  constante.  V oor  radar 
beelden  wordt  vaak  een  volgende  template  gebmikt  [Novak  1991,  1992]: 


M  1 1 1 1 1 1 1 1 1 1 1  r 

- 

- 

■ -JCf  □  ,  target 

M  1  1  1  1  1  i  1  1  1  1  IT 

clutter 


Figuur  4.4:  CFAR  template  met  target. 

De  betrouwbaarheid  van  het  gemiddelde  ^  en  de  standaard  deviatie  tJc  wordt 
bepaald  door  de  omvang  van  de  template  (het  aantal  achtergrond  samples).  Hoe 
hoger  het  aantal  samples  des  te  lager  het  detectie  verlies.  De  omvang  van  de  tem¬ 
plate  is  afhankelijk  van  de  omvang  van  de  doelen,  de  omgeving  en  de  pixelafme- 
tingen.  Wordt  een  template  te  klein  ten  opzichte  van  het  doel,  dan  bevinden  zich  te 
veel  doelpixels  tussen  de  clutterpixels  waardoor  het  gemiddelde  en  de  standaard 
deviatie  minder  betrouwbaar  zijn.  Wordt  een  template  te  groot  dan  bestaat  de 
mogelijkheid  dat  een  doel  niet  meer  wordt  vergeleken  met  zijn  directe  achtergrond, 
maar  met  een  andere,  verder  weg  gelegen  achtergrond.  In  bijvoorbeeld  [Novak 
1991]  is  gekozen  voor  een  template  van  41  x  41  pixels  (totaal  160)  corresponde- 
rend  met  een  afmeting  van  12.5  m  x  12.5  m  (pixelafmeting  0.30  m  x  0.30  m)  ten 
behoeve  van  het  detecteren  van  stationaire  militaire  doelen.  De  CFAR  constante  is 
afhankelijk  van  de  gewenste  kans  op  een  false-alarm  veroorzaakt  door  de  ruis  in 
het  radar  beeld.  De  relatie  is: 
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jpix\l  =  l)dx  (21) 

Hier  is  de  kans  op  een  false-alarm,  x  de  intensiteit  of  het  vermogen  van  een 
sample,  I  de  onderliggende  intensiteit  zonder  speckle  (het  gemiddelde  van  x),  en 
P(;cl/ )  de  gamma  verdeling  van  x  gegeven  I  [de  Vries  1995].  De  relatie  tussen  de 
standaard  deviatie  en  het  gemiddelde  van  x  wordt  gegeven  door  vergelijking 
(5).  Dit  geldt  zowel  voor  de  oorspronkelijke  kanalen  HH,  HV  en  W,  als  voor  de 
resultaten  van  de  polarimetrische  filters  gegeven  in  paragraaf  4.1. 

Uit  vergelijking  (20)  blijkt  dat  de  CFAR  detector  gebaat  is  bij  een  zo  hoog  moge- 
lijke  target-to-clutter  ratio  (t/c  ratio:  Xt  -  fA)  en  een  zo  klein  mogelijke  standaard 
deviatie  (CTc)-  Anders  gezegd;  een  CFAR  detector  werkt  het  best  wanneer  voor  een 
doel  geldt  dat  de  t/c  ratio  gedeeld  door  de  standaard  deviatie  (deflection  ratio  ) 
maximaal  is.  De  in  paragraaf  4.1  beschreven  filters  als  het  PWF  reduceren  alleen 
de  standaard  deviatie  van  de  achtergrond.  Er  bestaan  ook  filters  die  ontworpen  zijn 
om  de  t/c  ratio  te  vergroten  [Novak  1989,  Verbout  1992].  Het  PWF  levert  echter 
een  hogere  deflection  ratio  op  dan  al  deze  filters  [Novak  1992,  Verbout  1992]  en  is 
eenvoudiger  van  opzet.  Dit  is  ook  de  reden  waarom  deze  filters  niet  zijn  uitgewerkt 
in  paragraaf  4.1. 


Order-statistics  CFAR  detector 

Bij  order-statistics  zijn  de  samples  van  een  verdeling  gerangschikt  naar  de  grootte 
van  de  waarde  [Rohling  1983].  Zo  is  bijvoorbeeld  X30  de  waarde  die  zich  op  30% 
van  het  totale  aantal  bevindt.  Een  order-statistics  CFAR  detector  wordt  beschreven 
door  [Novak  1991]: 


- ->K, 


CFAR 


(22) 


Hier  is  x,  de  sample  die  getest  wordt,  jc„  de  mediaan  (een  benadering  voor  het 
gemiddelde,  ook  wel  aangeduidt  met  X50)  en  Xp  -  een  benadering  van  de  stan¬ 
daard  deviatie  (bijvoorbeeld  xeo  -  .^40)-  Het  voordeel  van  order-statistics  CFAR 
detectoren  is  dat  zij  robuuster  zijn  in  het  geval  dat  zich  doelpixels  tussen  de  clut- 
terpixels  bevinden.  Dit  kan  voorkomen  wanneer  de  template  noodgedwongen  klein 
is  of  wanneer  zich  meerdere  doelen  naast  elkaar  bevinden.  De  order-statistics 
CFAR  detector  kent  twee  varianten:  (a)  de  hollow-stencil  variant  waarbij  de  tem¬ 
plate  er  uit  ziet  als  in  figuur  4.4  en  (b)  de  solid-stencil  variant  waarbij  alle  pixels 
die  zich  bevinden  tussen  de  rand  en  pixel  x,  worden  meegenomen  in  de  berekening 
van  het  gemiddelde  en  de  standaard  deviatie. 


Weibull  detector 

De  pixel  amplitudes  van  een  homogeen  gebied  binnen  een  SAR  beeld  zijn 
Rayleigh  verdeeld  [Ulaby  1982].  Bevindt  zich  in  dit  gebied  een  doel  waardoor  het 
niet  meer  homogeen  is,  dan  verandert  die  verdeling  in  een  Weibull  verdeling 
[Carlotto  1996]: 


(23) 


a 

p{x)  =  — 
m 


A  v^ 


a-l 


exp 


\mj 


Hierin  is  ;c  de  amplitude  van  een  sample,  m  een  schaal  parameter  en  a  de  Weibull 
parameter.  Opgemerkt  dient  te  worden  dat  in  (21)  x  de  intensiteit  (het  vermogen) 
voorstelt,  het  kwadraat  van  de  amplitude.  De  pixel  intensiteiten  zijn  overigens  niet 
Rayleigh  verdeeld.  Vergelijking  (23)  geldt  uitsluitend  voor  1-look  beelden.  De 
Rayleigh  verdeling  is  die  van  (23)  voor  cc  =2.  Figuur  4.5  toont  de  Weibull  verde- 
ling  voor  m  =  2  en  verschillende  waarden  van  a. 


X 


Figuur  4.5:  De  Weibull  verdeling  voor  m  =  2en  verschillende  waarden  van  a. 

Figuur  4.5  laat  zien  dat  wanneer  a  =  1,  de  Weibull  verdeling  exponentieel  is.  Wat 
opvalt  is  dat  naarmate  a  groter  wordt,  de  verdeling  symmetrischer  wordt  en  het 
maximum  verschuift  naar  m.  Voor  de  mediaan  geldt: 

X  =ln“(2)m  (24) 

Het  gemiddelde  en  de  variantie  (het  kwadraat  van  de  standaard  deviatie)  zijn 
moeilijker  af  te  leiden.  Tabel  4.3  toont  de  gemiddelden  en  varianties  voor  enkele 
waarden  van  a  (analytische  uitdrukkingen  worden  gegeven  in  [Levanon  1988]). 
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Tabel  4.3:  Het  gemiddelde  en  de  variantie  van  de  Weibull  verdeling  voor  enkele  waar- 
den  van  a. 


a 

gemiddelde 

variantie 

0.25 

24m 

39744 

0.5 

2m 

20 

1 

m 

2 

{i-fK 

In  [Carlotto  1996]  wordt  voor  een  SAR  beeld  van  512  x  512  pixels  a  uitgerekend 
voor  gebieden  of  subbeelden  van  8  x  8  pixels.  Dit  wordt  gedaan  door  het  cumula- 
tieve  histogram  van  de  subbeelden  te  vergelijken  met  vooraf  uitgerekende  cumula- 
tieve  histogrammen  (voor  !<«  <4,  in  32  stappen).  Het  uitgerekende  cumulatieve 
histogram  dat  het  minste  verschilt  van  dat  van  het  subbeeld  geeft  a.  Het  blijkt  dat 
een  histogram  van  de  uitgerekendecr ’s  bimodaal  is,  dat  wil  zeggen,  twee  maxima 
vertoont,  gescheiden  door  een  minimum  (bij  CX  =2.7).  Het  deel  rechts  van  het 
minimum  wordt  geassocieerd  met  natuurlijk  terrein,  het  deel  links  hiervan  met 
niet-natuurlijk  (man-made)  terrein.  Bebouwde  gebieden  en  een  landingsbaan 
worden  op  deze  manier  gedetecteerd. 

De  bovenstaande  manier  van  bepaling  van  (X  is  wat  omslachtig.  Uit  figuur  4.5 
blijkt  dat  de  variantie  van  de  Weibull  verdeling  afneemt  met  toenemende  a  (de 
verdeling  wordt  “spitser”).  Het  is  dus  ook  mogelijk  a  te  schatten  door  de  variantie 
te  berekenen.  Figuur  4.6  laat  zien  welk  verband  er  is  tussen  a  en  de  variantie  voor 
m=l. 


Figuur  4.6:  Het  verband  tussen  a  en  de  variantie  van  de  Weibull  verdeling. 

De  curve  is  berekend  door  numerieke  integratie  voor  Of  =  0.50,  0.55,  0.60,  0.65, ... 
5.00).  Een  moeilijkheid  is  dat  er  twee  onbekende  parameters  zijn:  a  en  m .  Met 
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behulp  van  vergelijking  (24),  het  feit  dat  de  variantie  recht  evenredig  is  met 
(tabel  4.3)  en  een  eenvoudig  iteratief  algoritme  zijn  ze  beide  te  bepalen.  In 
[Carlotto  1996]  wordt  dit  probleem  vermeden  door  onterecht  aan  te  nemen  dat  m 
de  mediaan  is.  Deze  benadering  is  vooral  voor  kleine  a  onjuist,  zie  vergelijking 
(24).  Een  derde  manier  omcc  te  berekenen  wordt  gegeven  in  [Sekine  1992].  Hierin 
wordt  het  probleem  gereduceerd  tot  een  kleinste  kwadraten  fit  van  een  rechte  lijn 
aan  de  pixelwaarden  van  een  subbeeld.  De  drie  methoden  (“Sekine”,  “dit  rapport” 
en  “Carlotto”)  zijn  met  gesimuleerde  Weibull  verdeelde  data  getest.  Figuur  4.7 
toont  de  resultaten  voor  m=10  ena  =  0.50,  0.51,  0.52,  ...,  4.00.  De  drie  curven 
zijn  het  verschil  tussen  de  berekende  en  de  ingestelde  waarden  (er  is  een  offset  van 
2, 4  en  6  toegepast).  Een  maat  voor  de  nauwkeurigheid  van  een  methode  is  de  som 
van  de  absolute  verschillen  uit  de  figuur.  De  genormeerde  sommen  zijn  1.2 
(Sekine),  1.0  (dit  rapport)  en  2.1  (Carlotto).  De  methode  voorgesteld  in  dit  rapport 
is  dus  het  nauwkeurigst.  Bovendien  is  het  de  snelste  methode:  de  genormeerde 
rekentijden  zijn  5.3  (Sekine),  1.0  (dit  rapport)  en  14  (Carlotto). 


Figuur  4. 7:  Vergelijking  van  de  drie  methoden  om  a  te  schatten. 

Resultaten 

Bij  de  implementatie  van  de  CFAR  detectoren  is  gekozen  voor  een  template  van 
21x21  pixels,  ongeveer  32  m  x  32  m.  Deze  template  bleek  het  beste  compromis 
tussen  de  omvang  van  de  doelen,  de  omgeving  en  de  pixelafmeting.  Aangenomen 
werd  dat  een  totaal  aantal  pixels  van  80  een  voldoende  betrouwbaar  resultaat  zou 
opleveren  voor  het  gemiddelde  en  de  standaard  deviatie.  Dit  aantal  geldt  voor  de 
hollow-stencil  detectoren  waaronder  de  klassieke  CFAR  detector.  Een  solid-stencil 
van  21  X  21  pixels  bevat  totaal  440  clutterpixels.  Bij  de  order-statistics  detectoren 
is  nagegaan  welke  order-statistics  de  beste  benaderingen  gaven  voor  het  gemiddel¬ 
de  en  de  standaard  deviatie.  Voor  de  hollow-stencil  variant  bleek  het  beste  criteri- 
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^50 
^85  ~  ^15 


>K. 


CFAR 


Voor  de  solid-stencil  variant  bleek  dit: 


(25) 


>K, 


CFAR 


(26) 


De  afstand  in  de  noemer  van  de  solid-stencil  variant  is  kleiner  dan  bij  de  hollow- 
stencil  variant.  Dit  laat  zich  verklaren  door  het  feit  dat  de  clutterpixels  van  de 
solid-stencil  variant  meer  doelpixels  bevatten  bij  de  aanwezigheid  van  een  doel  in 
het  centrum.  Voor  alle  CFAR  detectoren  is  gekozen  voor  Kcfar  =  9.5  omdat  deze 
visueel  het  beste  overzicht  gaf  van  de  doelen.  Met  behulp  van  (21)  en  de  gamma 
intensiteit  verdeling  [de  Vries  1995]  is  uitgerekend  wat  de  false-alarm  kansen  zijn 
voor  verschillende  number-of-looks: 


label  4.4:  De  kans  op  een  false-alarm  voor  verschillende  number-of-looks,  Kcfar  -  9.5. 


ENL 

Pfa 

percentage 

3 

10® 

0.0001  % 

8 

10'® 

0.000001  % 

20 

10-10 

0.00000001  % 

Opgemerkt  dient  te  worden  dat  de  waarden  van  tabel  4.4  zijn  berekend  op  basis 
van  de  speckle.  Het  werkelijke  aantal  false-alarms  is  hoger  omdat  er  ook  objecten 
worden  gedetecteerd  waarin  men  niet  is  geinteresseerd.  Uit  tabel  4.4  blijkt  dat  de 
kans  op  een  false-alarm  afneemt  naarmate  het  aantal  looks  toeneemt.  Vergelijken 
we  deze  resultaten  met  die  uit  tabel  4.2  dan  zien  we  bijvoorbeeld  dat  het  PWF  de 
kans  op  een  false-alarm  ongeveer  100  x  verkleint. 


Bijlage  C  toont  de  resultaten  van  de  klassieke  CFAR  detector  (20)  en  de  order- 
statistics  detectoren  (25)  en  (26),  alle  drie  met  Kcfar  =  9.5,  voor  alle  filters  be- 
sproken  in  paragraaf  4.1.  Hierbij  zijn  de  gedetecteerde  pixels  afgebeeld  op  de 
originele  maar  donkerder  gemaakte  beelden.  Ook  hieruit  blijkt  dat  de  kans  op  een 
false-alarm  afneemt  naarmate  het  aantal  looks  toeneemt.  Vergelijken  we  beide 
hollow-stencil  CFAR  detectoren  dan  zien  we  dat  de  order-statistics  detector  min¬ 
der  false-alarms  genereert  dan  de  klassieke  detector.  De  order-statistics  detector  is 
dan  ook  robuuster,  maar  de  oorzaak  kan  ook  liggen  in  de  keuze  van  benadering  van 
de  standaard  deviatie.  Hierdoor  kan  de  kans  op  een  false-alarm  voor  deze  detector 
lager  zijn.  Het  verschil  tussen  beide  order-statistics  CFAR  detectoren  is  klein.  De 
solid-stencil  variant  lijkt  iets  beter  te  werken  dan  de  hollow-stencil  variant.  Een 
oorzaak  hiervan  kan  zijn  het  verschil  in  het  aantal  achtergrond  samples 
(respectievelijk  440  en  80)  en  daarmee  de  betrouwbaarheid  van  het  gemiddelde  en 
de  standaard  deviatie.  Een  waarschijnlijker  oorzaak  is  het  verschil  in  afstand  in  de 
noemer  van  vergelijking  (25)  en  (26)  als  benadering  voor  de  standaard  deviatie. 
Het  sorteren  van  440  pixels  (order-statistics)  duurt  wel  langer  dan  van  80  pixels. 
Op  basis  van  de  resultaten  gaat  de  voorkeur  uit  naar  de  order-statistics  detectoren. 
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In  bijlage  C  zien  we  ook  dat  het  Wavelet  Shrinkage  filter  het  beeld  enigszins  heeft 
aangetast  zodat  een  enkel  doel  niet  meet  gedetecteerd  wordt.  Het  meest  optimale 
filter  voor  detectie  lijkt  het  PWF  zijn.  Het  DIWF  is  een  goede  tweede  maar  geldt 
specifiek  voor  het  dihedral  beeld.  Om  na  te  gaan  wat  een  CFAR  detector  precies 
detecteert  is  er  ook  op  pixel  niveau  naar  deze  beelden  gekeken.  Figuur  4.8  toont 
een  uitvergroting  van  het  resultaat  van  de  hollow-stencil  order-statistics  CFAR 
detector  op  de  PWF  beelden  van  de  PHARUS  opnamen  180X  en  270X.  In  de 
eerste  opname  worden  vier  van  de  vijf  doelen  gedetecteerd,  in  de  tweede  alle  vijf 
de  doelen.  Het  aantal  false-alarms  in  beide  is  ongeveer  gelijk.  De  meeste  ervan 
blijken  te  worden  veroorzaakt  door  bomen,  zie  Bijlage  A.  Op  basis  van  figuur  4.8 
kan  niet  direct  worden  bepaald  welke  objecten  (een  groep  ogenschijnlijk  bij  elkaar 
horende  pixels)  een  voertuig  zijn  en  welke  niet.  De  pixels  zijn  echter  goed  gegroe- 
peerd  en  er  kan  een  uitspraak  worden  gedaan  over  de  afmeting  en  soms  ook  over 
de  richting  van  een  gedetecteerd  object. 

Voor  de  Weibull  detector  zijn  de  PHARUS  opnamen  opnieuw  verwerkt,  zonder 
multi-looking,  omdat  deze  verdeling  alleen  geldt  voor  1-look  data.  De  a  -waarden 
werden  bepaald  in  gebieden  van  5  x  5  pixels,  met  behulp  van  de  eerder  beschreven 
optimale  methode  (“dit  rapport”).  Figuur  4.9  toont  het  histogram  van  a  voor  het 
HH-beeld  van  opname  180X.  De  verdeling  lijkt  een  licht  scheve  normale  verde¬ 
ling,  met  het  maximum  bij  a  =2.  Dit  laatste  betekent  dat  de  pixelwaarden  Rayleigh 
verdeeld  zijn,  zoals  valt  te  verwachten  voor  een  1-look  HH  beeld.  Er  is  geen  sprake 
van  een  bimodale  verdeling  zoals  in  [Carlotto  1996],  waarschijnlijk  doordat  onze 
beelden  maar  zeer  weinig  doelen  bevatten.  Bijlage  C  bevat  de  resultaten  van  de 
Weibull  detector  voor  HH,  HV  en  W.  De  gedetecteerde  man-made  gebieden  zijn 
op  eenzelfde  manier  weergegeven  als  bij  de  CFAR  detectoren.  Hierbij  zijn  gebie¬ 
den  van  5  X  5  pixels  met  een  a  lager  dan  0.75-1.05  als  doelen  aangemerkt.  De 
drempel  werd  zo  gekozen  dat  er  ongeveer  10  doelen  gedetecteerd  werden.  Te  zien 
valt  dat  geen  enkel  doel  door  de  Weibull  detector  wordt  gedetecteerd.  Er  worden 
alleen  maar  false-alarms  waargenomen.  De  Weibull  detector  is  dus  ongeschikt 
voor  het  detecteren  van  doelen  met  PHARUS  data.  Wat  opvalt  aan  de  1-look 
beelden  is  dat  er  meerdere  puntdoelen  te  onderscheiden  zijn  binnen  de  grenzen  van 
de  voertuigen,  zie  ook  figuur  4.8.  Omdat  op  deze  data  geen  multi-looking  is  toege- 
past  is  de  azimuth  resolutie  1.2  m  gebleven,  zie  ook  hoofdstuk  3. 
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1=  GECKO 
2= TEAR 

3  =  ASTROSS 

4  =  MTLB 

5  =  ZSU  23-4 

Figuur  4.8:  De  resultaten  van  de  hollow-stencil  order-statistics  CFAR  detector  op  PWF 
data. 
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Figuur  4.9:  Histogram  van  a  van  opname  180X  HH. 


4.3  Herkenning 

Herkenning  is  het  maken  van  onderscheid  tussen  de  gedetecteerde  doelen.  Daarbij 
gaat  het  om  het  onderscheid  tussen  stationaire  militaire  voertuigen  en  false-alarms 
en  tussen  de  voertuigen  onderling.  Belangrijk  hierbij  zijn  de  kenmerken  of  featu¬ 
res,  op  basis  waarvan  een  gedetecteerd  doel  geclassificeerd  kan  worden.  In  deze 
paragraaf  worden  twee  methoden  uitgewerkt.  Hiertoe  is  op  basis  van  figuur  4.8  een 
verzameling  bepaald  van  bekende  doelen  en  false-alarms  bepaald  die  zoveel  mo- 
gelijk  in  beide  opnamen  (180X  en  270X)  voorkwamen,  zie  figuur  4.10. 

Omdat  de  pixels  per  gedetecteerd  object  goed  gegroepeerd  zijn,  kan  een  uitspraak 
worden  gedaan  over  de  afmetingen  en  soms  ook  de  richting  van  de  objecten.  Zij 
zijn  gegeven  in  tabel  4.5.  De  lengte  is  bepaald  door  het  aantal  pixels  op  te  meten  in 
de  ogenschijnlijke  orientatie  richting.  De  orientatie  is  ingedeeld  in  vier  klassen:  - 
45°,  0°,  45°  en  90°.  Er  wordt  dus  geen  onderscheid  gemaakt  tussen  de  voor-  en 
achterkant  van  een  voertuig.  Is  een  waarde  niet  gegeven  dan  is  het  object  niet 
gedetecteerd  (object  4  opname  180X)  of  te  klein  om  een  richting  te  kunnen  bepa- 
len.  De  identiteit  van  object  6  t/m  10  is  vastgesteld  met  behulp  van  de  ground-truth 
data,  zie  Bijlage  A.  Object  1 1  kwam  hier  niet  in  voor  maar  bleek  een  mobiel  ver- 
hoogd  platform  voor  het  maken  van  opnamen  met  radiometers  [Bartlett  1997]. 
Wordt  alleen  de  lengte  in  ogenschouw  genomen  dan  kunnen  we  alle  objecten  van 
minimaal  7  m  lang  beschouwen  als  een  voertuig.  Dat  dit  ook  mis  kan  gaan  bewijst 
object  5  in  opname  270X.  Kleinere  voertuigen  zullen  over  het  hoofd  worden 
gezien.  Bovendien  verschillen  de  lengtes  van  object  4  en  5  aanzienlijk  tussen  beide 
opnamen.  De  oorzaak  hiervan  is  de  aard  van  het  voertuig  en  de  belichtingsrichting 
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(loodrecht  op  de  vliegrichting)  hiervan;  enkele  graden  verschil  kan  al  een  volledig 
andere  signatuur  opleveren  [de  Vries  1998].  Een  verschil  in  orientatie  van  een 
object  tussen  beide  opnamen  kan  ook  een  criterium  zijn  om  een  object  uit  te  slui- 
ten  van  de  verzameling  voertuigen.  Dit  gaat  echter  wederom  mis  voor  object  5  in 
opname  270X.  Het  is  begrijpelijk  dat  met  de  minimale  informatie  die  het  CFAR 
beeld  oplevert  geen  betrouwbare  classificatie  gedaan  kan  worden.  Er  is  meer 
informatie  nodig. 

Een  tweede  methode  is  herkenning  op  basis  van  de  covariantie  matrix.  Het  PWF 
CFAR  beeld  van  figuur  4.8  of  4.10  wordt  hierbij  gebruikt  als  “cue”  voor  de  ge- 
middelde  covariantie  matrix  van  een  object.  Op  basis  van  deze  matrices  kan  onder- 
scheid  worden  gemaakt  tussen  dihedral/trihedral  stracturen  en  natuurlijke  objec- 
ten.  Voertuigen  en  andere  kunstmatige  (man-made)  objecten  bezitten  vaak  meerde- 
re  rechte  hoeken  die  een  dihedral  of  trihedral  kunnen  vormen.  label  4.6  toont  de 
gemiddelde  covariantie  matrices  van  dezelfde  objecten.  Zij  zijn  verkregen  door 
rond  het  ogenschijnlijke  middelpunt  van  een  object  3x3  covariantie  matrices  te 
middelen.  Dit  is  nodig  omdat  vooraf  geen  onderscheid  gemaakt  kan  worden;  de 
covariantie  matrix  van  een  natuurlijk  object  is  meestal  een  stochastische  variabele 
met  een  gemiddelde  en  standaard  deviatie.  Opvallend  voor  dihedrals  en  trihedrals 
is  dat  de  HV  component  (e)  0  of  in  ieder  geval  erg  klein  moet  zijn.  Uit  tabel  4.6 
blijkt  dat  dit  ook  nogal  eens  voorkomt  bij  bomen.  Een  ander  kenmerk  voor  dihe¬ 
drals  en  trihedrals  is  de  fase  van  de  correlatie  coefficient  (p ).  Voor  een  goed 
gecalibreerd  beeld  horen  zij  zeer  uitgesproken  te  zijn.  Uit  tabel  4.6  blijkt  dit  ook 
weer  bij  bomen  voor  te  komen.  De  uitgesprokenheid  correspondeert  overigens  wel 
met  de  compositie  van  het  DIWF  en  TRIWF  in  figuur  4.3.  Vooral  in  opname  270X 
is  duidelijk  te  zien  dat  object  1,  6,  en  11  een  dihedral  reflectie  aangeven  (duidelijk 
rode  kleur  die  correspondeert  met  DIWF).  De  resterende  covariantie  matrix  para¬ 
meters  doen  ook  weinig  verschil  vermoeden  tussen  de  voertuigen  en  bomen.  Een 
oorzaak  van  het  kleine  verschil  tussen  de  covariantie  matrices  kan  zijn  dat  het 
aantal  matrices  dat  wordt  gemiddeld  (3  x  3)  niet  voldoende  is  voor  een  betrouw- 
baar  gemiddelde.  Het  verkrijgen  van  meer  samples  en  een  betere  resolutie  is  echter 
niet  mogelijk  zonder  de  radar  aan  te  passen. 

Uit  de  twee  behandelde  voorbeelden  blijkt  in  ieder  geval  dat  classificatie  van  de 
doelen  moeilijk  is  bij  de  huidige  resolutie.  Door  gebrek  aan  tijd  binnen  dit  project 
is  niet  verder  gekeken  naar  andere  methoden.  Dit  is  wel  gebeurd  in  het  vervolg  op 
dit  onderzoek  [de  Vries  1998].  Hierbij  is  ook  gekeken  naar  data  met  een  hogere 
resolutie. 
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Figuur  4.10:  De  verzameling  gedetecteerde  objecten  (bekende  doelen  en  false-alarms)  die 
is  gebruikt  voor  het  testen  van  enkele  herkenningsmethoden.  De  nummers 
corresponderen  met  label  4.5.  Het  uitgangspunt  zijn  de  resultaten  van  de 
hollow-stencil  order-statistics  CFAR  detector  op  PWF  data  ( figuur  4.8). 
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Tabel  4.5:  De  berekende  lengtes  en  richtingen  van  de  gedetecteerde  objecten  1  t/m  II. 


target 

identity 

180X 

length  direction 

270X 

length  direction 

1 

GECKO 

10  m 

45° 

12  m 

45° 

2 

TEAR 

12  m 

-45° 

12  m 

-45° 

3 

ASTROSS 

8  m 

-45° 

13  m 

-45° 

4 

MTLB 

- 

- 

8  m 

90° 

5 

ZSU  23-4 

7m 

90° 

3  m 

0° 

6 

tree 

6  m 

45° 

6  m 

O 

1 

7 

tree 

6  m 

45° 

6  m 

90° 

8 

tree 

2  m 

- 

4  m 

-45° 

9 

tree 

5  m 

90° 

2  m 

- 

10 

tree 

5  m 

-45° 

5  m 

0° 

11 

platform 

11m 

-45° 

9m 

-45° 

Tabel  4.6:  De  3x3  gemiddelde  covariantie  matrices  van  de  gedetecteerde  objecten  1 
t/m  1 1.  Opgenomen  zijn  ook  de  matrices  van  een  dihedral  (a=0° )  en  een 
trihedral  [Freeman  1992]. 


target 

identity 

e 

r 

Ipl 

180X 

e 

y 

Ipl 

270X 

dihedral 

0 

1 

1 

180° 

0 

1 

1 

180° 

trihedral 

0 

1 

1 

0° 

0 

1 

1 

0° 

1 

GECKO 

0.27 

1.20 

0.35 

-114° 

0.13 

1.47 

0.84 

158° 

2 

TEAR 

0.02 

0.53 

0.91 

-27° 

0.14 

1.13 

0.58 

-7° 

3 

ASTROSS 

0.04 

0.54 

0.68 

-56° 

0.12 

1.77 

0.61 

66° 

4 

MTEB 

0.09 

1.31 

0.65 

-22° 

0.26 

5.59 

0.72 

-1° 

5 

ZSU  23-4 

0.03 

0.16 

0.66 

139° 

0.16 

0.65 

0.76 

90° 

6 

tree 

0.23 

0.49 

0.71 

-33° 

0.67 

1.90 

0.52 

129° 

7 

tree 

0.22 

1.24 

0.72 

-65° 

0.20 

0.81 

0.64 

-61° 

8 

tree 

0.17 

0.58 

0.39 

-17° 

0.05 

0.97 

0.79 

-56° 

9 

tree 

0.04 

0.33 

0.69 

31° 

0.86 

1.26 

0.77 

-3° 

10 

tree 

0.34 

1.41 

0.35 

1 

o 

o 

0.07 

0.79 

0.80 

-3° 

11 

platform 

0.31 

1.48 

0.81 

52° 

0.09 

0.34 

0.74 

-142° 
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5.  Condusies  en  aanbevelingen 


Binnen  dit  onderzoek  is  gekeken  naar  het  opsporen  van  stationaire  militaire  voer- 
tuigen  met  behulp  van  PHARUS,  een  polarimetrische  SAR.  Doelopsporing  bestaat 
uit  drie  stappen:  pre-processing,  detectie  en  herkenning.  In  de  pre-processing 
wordt  de  data  zoveel  mogelijk  ontdaan  van  speckle  mis. 

Diverse  polarimetrische  speckle  reductie  algoritmen  zijn  onderzocht.  Het  beste 
bleek  het  PWF  (Polarimetrisch  Whitening  Filter).  Het  DIWF  (Dihedral  Whitening 
Filter)  bleek  een  goede  tweede  met  betrekking  tot  de  detectie  van  dihedral  stmctu- 
ren.  Het  Span  en  OIS  (Optimal  Intensity  Sum)  presteerden  minder  maar  bleken  een 
goed  altematief  door  hun  eenvoud.  Het  enige  niet  polarimetrische  filter  dat  werd 
onderzocht,  het  Wavelet  Shrinkage  filter,  zorgde  voor  een  kleine  verbetering 
bovenop  het  PWF  maar  bleek  de  resolutie  aan  te  tasten. 

Opvallend  goed  waren  de  resultaten  van  de  CFAR  detectoren  op  PHARUS  data; 
alle  doelen  (op  een  doel  in  een  beeld  na)  werden  gedetecteerd  en  het  aantal  false- 
alarms  bleef  beperkt.  De  beste  detector  bleek  de  order-statistics  CFAR  detector. 

De  voorkeur  gaat  hierbij  uit  naar  de  “hollow  stencil”  variant  omdat  deze  sneller  is. 
De  Weibull  detector  miste  alle  doelen.  Het  maken  van  onderscheid  tussen  de 
militaire  voertuigen  onderling  en  tussen  false-alarms,  bleek  moeilijk.  Op  basis  van 
het  CFAR  beeld  kon  een  uitspraak  worden  gedaan  over  de  afmetingen  en  de  ori- 
entatie  van  een  gedetecteerd  object.  Volledige  classificatie  bleek  echter  niet  mo¬ 
gelijk,  ook  niet  op  basis  van  de  polarimetrische  eigenschappen  (covariantie  ma¬ 
trix).  De  belangrijkste  oorzaak  hiervan  is  de  resolutie  (circa  3.5  m  x  5.0  m). 

Het  detecteren  van  doelen  met  PHARUS  blijkt  dus  succesvol.  De  polarimetrische 
mode  draagt  vooral  bij  in  de  reductie  van  speckle.  Herkenning  vergt  echter  een 
hogere  resolutie  en  nader  onderzoek  naar  classificatie  methoden.  Een  deel  van  dit 
onderzoek  is  reeds  uitgevoerd  [de  Vries  1998],  waarbij  is  gekeken  naar  polarime¬ 
trisch  gefilterde  data  met  een  hogere  resolutie  dan  PHARUS  op  dit  moment  biedt. 
Aanbevolen  wordt  ook  te  kijken  naar  de  polarimetrische  eigenschappen  van  de 
hoge  resolutie  SAR  data  die  is  opgenomen  tijdens  dezelfde  proefnemingen 
[Bartlett  1997].  Deze  data  is  reeds  beschikbaar. 

Doelopsporing  kan  bovendien  het  beste  worden  uitgevoerd  met  meerdere  sensoren 
waarbij  PHARUS  aangeeft  waar  zich  mogelijke  doelen  bevinden  (cueing).  Met 
andere  sensoren  kan  dan  worden  onderzocht  wat  de  aard  van  deze  doelen  is.  Voor- 
beelden  van  goede  aanvullende  sensoren  zijn  thermisch-infrarood  (TIR),  laser- 
radar  en  multispectrale  scanners.  Het  voordeel  van  SAR  en  laser-radar  is  dat  zij 
ook  in  staat  zijn  hoogte  of  diepte  weer  te  geven,  informatie  die  van  groot  belang 
kan  zijn  bij  herkenning. 
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Bijlage  A  Ground-truth 


Figuur  A.l:  Luchtfoto  van  het  gebied  waar  de  voertuigen  zijn  opgesteld  ( zander  voertui- 
gen)  van  november  1995,  hoogte  6280  m. 
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Figuur  A.2:  Tweefoto 's  van  het  gebied  waar  de  voertiiigen  zijn  opgesteld  ( zander  voer- 
tuigen),  genomen  vanidt  een  helikopter  in  de  richting  noord-oost  (boven)  en 
noord-west  (onder). 
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Figuur  A.4:  Panorama  foto’s  van  het  gebied  waar  de  voertuigen  zijn  opgesteld  (met 
voertuigen). 
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Figuur  A.5:  Kaart  van  het  gebied  waar  de  voertuigen  zijn  opgesteld  ( met  voertuigen). 
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